Non-linear theory of deformable superconductors 



00 

o 
o 



Pavel Lipavsky 1,2 , Klaus Morawetz 3,4 , Jan Kolacek 2 and Ernst Helmut Brandt 5 
1 Faculty of Mathematics and Physics, Charles University, Ke Karlovu 3, 12116 Prague 2, Czech Republic 
2 Institute of Physics, Academy of Sciences, Cukrovarnickd 10, 16253 Prague 6, Czech Republic 
3 Forschungszentrum Rossendorf, PF 51 01 19, 01314 Dresden, Germany 
4 Max- Planck- Institute for the Physics of Complex Systems, 
Noethnitzer Str. 38, 01181 Dresden, Germany and 
5 Max-Planck-Institute for Metals Research, D-70506 Stuttgart, Germany 

The interaction of the superconducting condensate with deformations of the crystal lattice is 
formulated assuming the electrostatic potential to be of Bernoulli type and the effect of strain on 
material parameters. In the isotropic approximation it is shown that within the Ginzburg-Landau 
theory both contributions can be recast into the local but non-linear interaction term of the free 
energy. 
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PACS numbers: 74.20.De, 74.25. Ld, 74.81.-g 



I. INTRODUCTION 

When cooled, metals reduce their volume. At the 
transition from normal to superconducting state, the co- 
efficient of the thermal expansion makes a jump. In 
most cases, the superconducting systems reduce their 
volume less than the normal ones. Consequently, inho- 
mogeneities of the superconducting phase cause stresses 
which are similar to stresses caused by the inhomo- 
geneities due to the temperature. 

Since the superconducting condensate affects the spe- 
cific volume, deformations of the crystal lattice also affect 
the condensate. Mechanisms of this interaction between 
the crystal lattice and the condensate can be outlined 
within the simplest two-fluid free energy of Gorter and 
Casimir, /gc = — jlT^uj — ^jT 2 ^/l — u>, where ui is the 
superconducting fraction, 7 is the linear coefficient of the 
specific heat per unit volume known as the Sommerfeld 
7, and T c is the critical temperature. 

Material parameters 7 and T c are not constants. They 
depend on the electron density n and the deformation 
of the crystal lattice which we describe by the lattice 
density riiat- Due to the dependencies 7(71, m a t) and 
T c (n, Uiat) the crystal lattice interacts with the super- 
conducting condensate. 

In literature on deformable superconductors one 
finds two models of the lattice-condensate (lc) in- 
teraction. First, there are various phenomenological 
theoriesii^i 4 .! 5 -^ which assume that the density of elec- 
trons exactly follows the density of the lattice, n = n\ at - 
The perturbation of material parameters due to the 
deformation, e.g., £7 = (d^/dn)Sn + (dj/dni a t)5n} a t, 
then can be expressed via the lattice density, 5"/ — 
[(dj/dn) + (d^f/dni a t)] Sni at . The strength of the lc- 
interaction thus depends on the sum of both den- 
sity derivatives, (d^y/dri) + {d^/dn\ at ) and (dT c /dn) + 
(dT c /dn lat ). 

Second, a model in which the lc-interaction is medi- 
ated by the electrostatic potential has been discussed.-^ 
Since the theory of the electrostatic potential has been 
developed under the approximation of a stiff lattice, the 



lc-interaction obtained within this model depends exclu- 
sively on the derivatives with respect to the electron den- 
sity, dj/dn and dT c /dn. 

The phenomenological approach is more general be- 
ing applicable to all materials while the electrostatic ap- 
proach is limited to cases in which the dependence on the 
electron density dominates. On the other hand, the elec- 
trostatic approach offers a natural picture of the surface, 
in particular, one can easily see that the electrostatic field 
of the surface dipole contributes to the forces deforming 
the crystal lattice.— Studies within the phenomenological 
approach have not noticed the surface tension. 

In this paper we derive a phenomenological theory 
which unifies both approaches. To this end it is necessary 
to take into account the charge of a deformed lattice in 
the electrostatic potential and to allow for lc-interaction 
which is not covered by the (mean) electrostatic poten- 
tial. To avoid lengthy formulae or non-transparent tensor 
notation with numerous indices, we restrict our attention 
to the interaction between the lattice compression and 
the condensate. We neglect the interaction between the 
condensate anisotropy and shear deformations of the lat- 
tice. Interaction of the superconducting condensate with 
a deformation which can be interpreted as a mutual dis- 
placement of sublattices has been discussed in Ref. [lfj |. 



A. Origin of two mechanisms 

There are various microscopic mechanisms due to 
which material parameters of the superconductor depend 
either on the density of electrons or on the deformation 
of the crystal lattice. Although our discussion is indepen- 
dent of actual microscopic mechanisms, we find it prof- 
itable to outline some of these possibilities so that the 
need to treat perturbations of the electron density and 
the lattice density independently becomes more appar- 
ent. 

The Sommerfeld 7 is proportional to the single-spin 
density of states N at the Fermi energy, 7 = |7r 2 fc|Ao. 
The density of states naturally depends on the value of 
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the Fermi energy, which itself depends on the electron 
density. In this way the Sommerfeld 7 depends on the 
electron density. 

The density of states also reflects the electron band 
structure. For example, the saddle points giving a high 
density of states are quite sensitive to atomic spacing. 
Moreover, in ionic crystals of high-T c superconductors 
the charge transfer between sublattices depends on the 
lattice deformation. The Sommerfeld 7 thus depends on 
the lattice deformation via mechanisms which are distinct 
from changes of the electron density. 

The critical temperature is an even more complex 
quantity. For simplicity we express it within the BCS 
approximation T c = 0.85 9dcxp(— 1/NqV). Apparently, 
T c depends on the electron density and the deformation 
via the density of states. Beside this, there are addi- 
tional contributions via the interaction potential V and 
the Debye temperature 0d- For example, a compression 
of the lattice increases the mass density which results in a 
slower velocity of sound. This reduces the Debye temper- 
ature Od leading in some superconductors to a decrease 
of T c under an applied pressure ^ 

The above mentioned mechanisms of the density de- 
pendence work in pure materials. Let us mention a 
mechanism specific for dirty superconductors. In metals 
doped by paramagnetic impurities the dominant pressure 
dependence of T c results from the electron density depen- 
dence of the magnetic scattering relaxation time^ 



B. Plan of the paper 

In section [II] we introduce the free energy which com- 
bines the condensation energy of Ginzburg and Landau 
(GL), the energy of electric and magnetic fields, and the 
deformation energy. In section HTT1 we present the set of 
equations derived from the Lagrange variational princi- 
ple. In sections IIVI and [V] we focus on the electrostatic 
potential and the strain, respectively. In section I VII we 
write down an effective free energy for deformable super- 
conductors, and section fVIII is the summary. 



II. FREE ENERGY 

We start from the free energy and employ the La- 
grange variational principle to derive all stability con- 
ditions. Following Ginzburg and Landau (GL), the free 
energy of the superconducting state 

/. = fn + a\yj\ 2 + lp + -L |(_ ift v - e*A) 

(1) 

is defined as a non-local bi-quadratic function of the GL 
function ip. The non-local term has a form of the kinetic 
energy with the Copper pair mass m* and charge e* = 2e. 



The GL free energy is added to the free energy of the 
normal state 

/n = /o + e<p(n - nut) - U |V^| 2 + ^- |V x A| 2 . (2) 

The normal free energy covers the magnetic energy (last 
term), the electrostatic energy (second and third terms) 
and the local free energy f . 

The local free energy /o is a function of the electron 
density n, deviations of atomic positions u and temper- 
ature T 

/o = / (n,u,T). (3) 

According to the theory of elasticity^ the free energy 
does not depend directly on the vector u but only on its 
derivatives expressed via the strain tensor 

1 / dui dv,4 \ , 
»« = 2fe + *fJ- < 4 » 

There is no term explicitly attributed to the interaction 
between the deformation u and the superconducting con- 
densate, however this interaction has a number of hidden 
contributions. First, it is mediated by the electrostatic 
force. Second, the GL parameters a, 0, m* depend on 
the density of electrons and on the deformation u. 

For simplicity we will assume that the GL parameters 
a, (3 and m* depend on the strain exclusively via the 
crystal density ni at . In other words we neglect effects of 
shear deformations which break the isotropy of the sys- 
tem. The reader interested in the anisotropic interaction 
is referred to papers by Miranovic et a& or Cano et a%. 

In general, the lattice charge density in deformed crys- 
tal is a non-trivial problem since ni at describes the ionic 
charge and the deformation can change ionicity. We will 
assume that the ionicity is constant so that the charge is 
given by the divergence of atomic shifts 

n\ax = n (1 - (V • u)) = n (1 - u n - u 22 - u 33 ) . (5) 

Now all components of the free energy and material 
dependencies are specified. It remains to derive the equa- 
tions for the individual fields. To this end we employ the 
Lagrange variational principle. 

III. LAGRANGE VARIATIONAL CONDITIONS 

The free energy / s depends on the following indepen- 
dent variable fields: the vector potential A, the complex 
GL function the electrostatic potential <p, the electron 
density n, and the vector of atomic shifts u. The corre- 
sponding variations are well established , 14 ! 15 ' 16 therefore 
we present the resulting equations without derivations. 
We note that the vector and scalar potentials are in the 
Coulomb gauge, (V • A) = 0. 

The A-variation yields the Ampere law 

V 2 A = -^o-^Re ${-ifiS7 -e*A)ip. (6) 
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The ^-variation results in the GL equation 

f-iW - e*A) — !— {-ihV ~ e*A) ip + aip + /3\ip\ 2 i) = 0. 
2m* 

(7) 

The less usual form of the kinetic energy is a hermitian 
operator also for an inhomogeneous mass to*. The de- 
pendencies of the material parameters a, and to* on 
the lattice deformation u and the electron density n are 
rather weak but essential for specific problems like the 
vortex pinning by the strain around dislocations or for 
the effect of the elastic energy on the arrangement of 
vortices. 

The (/^-variation recovers the Poisson equation 



eVV = e(n - n lat ). 



(8) 



In deriving this equation we have neglected the ionic con- 
tribution to the dielectric function e, i.e., terms propor- 
tional to dt/dn\ a t. 

The n-variation furnishes us with the electrostatic po- 
tential known as the Bernoulli potential 



etp = - — - — \ip\ 2 
dn dn 

+ ip~(-ihV-e*A) 



1 d In to' 
2m* dn 



-ihV - e*A)tp. 

(9) 



The density derivative of the local free energy /o is non- 
trivial only if the system is perturbed from the homoge- 
nous state. To make dfo/dn transparent, it is necessary 
to expand it in perturbations. This rearrangement is ac- 
complished in section HVl 

The u- variation gives the strain equation 



dfo da 2 1 9/3 4 



dn 



lal 



du.ij " J duij 1 r 1 2" J duij 1 r 1 ' " J " r du; tJ 

— \ 111 777* — 

+ Vjtp (-ihV - e*A) 4> (-ihV - e*A) ip. 



2m* dui 



(10) 



We use the Einstein summation rule for doubled indices, 
e.g., Tjhjm = Y^j=i r jh'jm- The strain equation (fTU)) in- 
cludes terms which are so far rather symbolic. In sec- 
tion El we express all of them in terms of elastic moduli 
and forces on the crystal lattice. 



IV. BERNOULLI POTENTIAL 

We start with a rearrangement of the Bernoulli poten- 
tial © • The first derivative of the local free energy with 
respect to the electron density is the Fermi energy 



dfo 
dn 



= E 



F ■ 



(11) 



The Fermi energy itself depends on the electron den- 
sity via the Fermi-Dirac statistics and the exchange- 
correlation potential. — Besides, it depends on the lattice 



deformation via the density of states. Setting the Fermi 
energy of unperturbed system to zero, to the linear order 
in perturbations it reads 



E F 



dE F r dE F 
-on + — — u 



dn 



du^ 



(12) 



The Fermi energy E F depends on the lattice deforma- 
tion via changes of the electron band structure. Within 
the isotropic approximation we assume that it is propor- 
tional to the perturbation of the lattice density 



dE F dE F dniax 



dui 



dni a t dui 



dE F 
dn 



n Sj 



lal 



(13) 



Using the approximation JT5|) in the relation (fTTj) we ob- 
tain [ua = tin + U 2 2 + u 33 ] 



dfo dE F dE F 

— = -7^— on - n Ui 

dn dn croiat 



(14) 



The first term represents the Thomas-Fermi screening, 
the second one results from the charge inhomogeneity of 
the deformed ionic lattice. 



A. Thomas-Fermi Screening 

Now we express the change of the Fermi energy in 
terms of the electrostatic potential ip. To this end we 
use the Poisson equation (JSJ) in the form 



eV 2 ip = e(Sn + n^un). 



(15) 



Substituting Sn from equation (|15[) in the Fermi energy 
(fT4")l we arrive at 



dn dn e 



dE F dE F 



dn dn\ a 



noun. (16) 



The first term on the right hand side can be expressed 
via the Thomas-Fermi screening length 



dE F e 



dn e 2 TF 
The Bernoulli potential © now reads 



(17) 



ef — A 2 rF V 2 e(/3 = 



+ V(-iW-e*A) 



dE F dE F 
H — 

dn dnut 

da 2 19/3 

~2d^ 
1 d In to 



UqUu 



2m* dn 



(-ihV -e*A)ijj. 



(18) 



The electrostatic potential ip resulting from equation 
(US]) has two characteristic components, the free and the 
enforced one. The free solution is non-zero only near 
the surface decaying into the bulk on the Thomas-Fermi 
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screening length Atf- This solution is determined by a 
surface condition. We note that the free solution plays an 
important role in the surface dipole.— Here we focus on 
the bulk properties, therefore we ignore the free solution. 

Second, there is an electrostatic potential enforced by 
inhomogeneities in the superconducting density \%(j\ 2 and 
the lattice deformation as given by the right hand side of 
equation (JT5J) . We keep the name Bernoulli potential for 
this component. 

Two simplifications of the Bernoulli potential are at 
hand. First, we can neglect A TF V 2 e<,9. This is because 
gradients of the GL function and the corresponding po- 
tential are on the scale of the GL coherence length or the 
London penetration which are both much larger than the 
Thomas- Fermi screening length Atf- Second, the loga- 
rithmic derivative of the Cooper pair mass m* is a small 
quantity and we can neglect its gradient. Therefore 



etp : 



dE F 
dn 
dlnm* 
dn 



■0 (-ifiV 



da u a l d(3, ,4 
- e*A) — (-ihV - e*A) ip. 



2m 



(19) 



Note that neglecting the term V 2 ip in the Poisson 
equation (|15[) implies the quasi-neutral approximation 
n = ^lat- In this sense we can work with the non-zero 
electrostatic potential (fP9"]l while using the local charge 
neutrality for perturbations of material parameters. 

The Bernoulli potential (fT9|) extends previous 
result s 16 i 18 having two additional contribution. First, the 
charge of the deformed ion lattice is represented by the 
term oc tiqUu = rto(V-u). Second, the effect of the charge 
perturbation on the Cooper pair mass m* is included. 



B. From non-local to non-linear corrections 

The GL equation (0 multiplied by the conjugate GL 
function tp 

i> {-ihV - e*A) — !— (-ifiV - e*A) ib = -a|V>| 2 -/?|V>| 4 , 
2m* 

(20) 

couples the non-local term on the left hand side with the 
non-linear one /3|?/>| 4 . This gives us the freedom to make 
the Bernoulli potential either a linear or a local function 
of the superconducting density \ip\ 2 - We prefer the local 
but non-linear form, 



( dE-p dE ¥ 



\ dn dn 



n Ui 



da dhum* . 



(21) 



Apparently, there are a number of possible additional 
rearrangements of the Bernoulli potential. Since we 



study the interaction between the superconducting con- 
densate and the lattice deformation mediated by the 
Bernoulli potential, the form (j2"Tj) is optimal as it is ex- 
pressed in terms of Uu and \ip\ 2 - 



V. STRAIN EQUATION 



The strain equation (|10[) is rather involved as it con- 
tains gradients of derivatives with respect to tensor com- 
ponents of the strain. The major simplification follows 
from the assumption that all material parameters related 
to the superconducting phase depend on the strain ex- 
clusively via the lattice density ni a t, i-e., 



dc 



dui 



da 
dny 



-n Sj 



(22) 



and similar for (3 and m* . Within this isotropic approxi- 
mation the strain equation can be rearranged in a manner 
which in many steps parallels the treatment of the Fermi 
energy in the previous section. 



A. Stress 

The stress tensor has a general form of 

Pji = ^ 3 iikUki- (23) 

The moduli matrix A has 81 elements, but only 27 of 
them are independent^ 

Now we express the moduli tensor A in terms of the 
free energy /. We start with the strain-derivative of the 
local free energy 

df d 2 f d 2 f 

-um + —on, (24) 



duij duijduki duijdn 

which we have expanded in perturbations. The second 
term of expansion ([24]) we can express with the help of 
the already specified strain derivative of the Fermi energy 



d 2 f dE ¥ 



du^ dn duij 



dE F 
dn u 



■n S 



2i- 



(25) 



Finally we use the Poisson equation (|15p to eliminate the 
perturbation of the electron density Sn from the stress 

dfo d 2 f dEp 2 

a = "5 Z Ukl + OjiUkk-^ n 0- ( 26 ) 

duij duijduui dnut 

We have neglected the term Sji-^^n ^-V 2 ip, because it 
is proportional to A TF V 2 </?. 

An additional contribution to the stress results from 
the Coulomb interaction of the ionic lattice with itself. 
To make it explicit, we have to rearrange the electrostatic 
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term of the strain equation ([10]) with the help of the 
Bernoulli potential (|21|) 



eip-^— = -0 Jl n eip 



dui 



,'dE F dE F 

.da d In to* > 
-ouno \ — + a- 



dn dn 



Ud ^- dl ^)WV- <->7, 



2 V an 



The first term represents the Coulomb interaction of the 
lattice with itself. Other terms represent the interaction 
of the lattice with the superconducting condensate. 

The stress tensor collects all contributions to the strain 
equation (|10[) which arc linear in the strain u. The moduli 
matrix thus reads 



A 



jilk 



d 2 fo , t s ( n dE F , dE F , 2 



du rj duki 



+ S Tl S lk 2—^- + 
1 on- 



lat 



9n 



"o- (28) 



Here the second term arises from the increase of the elec- 
tron liquid energy under a volume compression. 



B. Deforming force 



In terms of the stress (|23|) the strain equation ([T 
reads 



V jPji — Ei . 



(29) 



where 



^_ , - .'da <3mm* \ da \ . .., 



1 {8(3 d In m* 



V,(^ - I ^ + 2/3- ^ 



i a/3 



2 



r) 1 ri 777 _ j _ 

+Vj — V - e*A) — 7 V HW ~ e*A) V 



5ui 



2m* 



(30) 



is the force (per unit volume) deforming the crystal. We 
have neglected the gradient of 9 In to* 

In the isotropic approximations (|22[) the deforming 
force (|3T)|) simplifies to a gradient 



F, = -n ViC/ 
of the effective potential 



(31) 



U 



da da d In to* 9 In to* , . 

a— I \wY 



dn dn 



lat 



dn 



dn ia 



l(W + JL + 2/3^ + 2/3^-'] \^ 



2 V dn dn 



lat 



Sri 



(32) 



We have used (|20|) to replace the non-local term of (|30|) 
by the non-linear one. 

As one can see, the force is given by the bi-quadratic 
effective potential with two material parameters 



b = 



da 
dn lat 

dp 
dn lat 



Oa 
dn 
dp 
dn 



dim 



d In? 



2/3 



dn ut 
d In to* 



dn 



■2/3 



lat 



dn 

d In to* 
dn 



(33) 



In both terms the derivatives enter in the same way as if 
one takes the volume or density derivative assuming the 
strict local charge neutrality. 



C. Isotropic model 

The simplest and mostly used isotropic model uses 
only two elastic moduli. The bulk modulus K measures 
changes of the specific volume and shear modulus u is the 
only coefficient of all volume- keeping deformations. For 
the isotropic system the strain equation (|29p simplifies 



t. 



,13 



K + -u) V(V.u) -/iVxVxu 



(34) 



where the force acting on a unit volume of the lattice is 
given by the gradient as 



F = aS7\ijj\ 2 + ^V|Vj| 4 . 



(35) 



Together with ([33]) this is our final result for the strain 
equation. 



VI. EFFECTIVE FREE ENERGY 

For studies of the lattice deformations it is not neces- 
sary to evaluate the electrostatic potential. In this case 
one can use a simplified free energy 

f' s = a\^ + i/3H 4 + -L |(- ifi V - e *A)V| 2 
2 2m* 
1 2 1 

+ -z — | V x A| + -A ijkl UijU k i 
zu I 



au ti \ik\ 2 - -bu u \ip\ 4 . 



(36) 



This free energy depends on the vector potential A, the 
GL function ip, and the displacement u. All material 
parameters a, /3, to*, A, a and b are now constant in 
space and do not undergo variations. 

By the Lagrange variation of the free energy f' s with 
respect to the vector potential A one recovers the Am- 
pere law ([5]). The variation of f' s with respect to the 
displacement u yields the strain equation (f2"9"]l with the 
force 432). 
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The effective free energy / s ' is not exactly equivalent to 
the full free energy / s , however. By the variation of f 
with respect to the GL function ip one arrives at the GL 
equation 

-i- (-iftV - e*A) 2 i>+(a-a uu)i>+(P-bu u )\i>\ 2 ^ = 0. 
2m* 

(37) 

Unlike the full GL equation ([7|), here the strain effect on 
the Copper pair mass m* is absent. It is mimicked by 
the m*-part of the strain effect on the effective potential, 
see a and b as given by equations (|33j) . 

VII. SUMMARY 

Starting from the free energy of the GL type we have 
derived the force which deforms the crystal lattice in 
the presence of the inhomogeneous superconducting con- 
densate. Neglecting terms proportional to the square of 
the small Thomas-Fermi screening length, we have rear- 
ranged the deforming force into the gradient of the bi- 
quadratic function of the GL function. 

Although we took into account perturbations of the 



charge neutrality and included the electrostatic potential, 
our result has confirmed that the assumption of the strict 
local charge neutrality can be applied for the evaluation 
of the force deforming the lattice. 

Based on our results, we have proposed an effective 
free energy which is simpler in being independent of the 
electrostatic potential and the density of normal elect- 
rons. Moreover, all its field variables are explicit so that 
there are no hidden interaction mechanisms. In particu- 
lar, it has no strain effect on the Copper pair mass m*. 
Contributions of these eliminated variables and depen- 
dencies arc covered by the effective local but non-linear 
interaction. 
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